function [Q,sse]=drawVarianceFull(nu,Q_0,T_0)
% input:
%       - yt is K (or N)x T
%       - B is K (or N) X K*p (or N*q)
T_1= T_0+size(nu,2);


sse=nu*nu';
Q_1=Q_0 + sse;
Q = iwishrnd(Q_1,T_1);
Q = diag(diag(Q));



